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Vorticity evolution in a rigid pipe of circular cross-section 


By F. Lam 


Yet every exact solution of the equations of fluid mechanics can actually 
occur in Nature. The flows must not only obey the equations which are 
non-linear, but also evolve in the distinct dynamic states as observed. 


In this paper, we show that the spatio-temporal evolution of incompressible 
flows in a long circular pipe can be described by vorticity dynamics. The principal 
techniques to obtain solution are similar to those used for flows in As the 
consideration of the Navier-Stokes equations is given in a cylindrical co-ordinates, 
two aspects of complication arise. One is the interaction of the velocity components 
in the radial and azimuthal directions due to the fictitious centrifugal force in the 
equations of motion. The rate of the vorticity production at pipe wall depends on 
the initial data at entry and hence is unknown a priori; it must be determined 
as part of the solution. The vorticity solution obtained defines an intricate flow 
field of multitudinous degrees of freedom. As the Reynolds number increases, the 
analytical solution predicts vorticity-scale proliferations in succession. For sufficient 
large initial data, pipe flows are of turbulent nature. The solution of the governing 
equations is globally regular and does not bifurcate in space or in time. It is asserted 
that laminar-turbulent transition is a dynamic process inbred in the non-linearity. 
The presence of exogenous disturbances, due to imperfect test environments or 
purpose-made artificial forcing, distorts the course of the intrinsic transition. The 
flow structures observed by Reynolds (1883) and others can be synthesised and 
elucidated in light of the current theory. 
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1. Introduction 

In the Eulerian description of the motion of an incompressible, homogeneous New¬ 
tonian fluid, the momentum and the continuity equations for fluid dynamics are 

9u/9t-I-(u.V)u = z/Au — V.u = 0, (1.1) 

where the velocity vector u = u(x, t) is the velocity, the scalar quantity p = p(x, t) 
is the pressure, the space variable is denoted by x=(a;, y, z), and A is the Laplacian. 
The density and the viscosity of the fluid are denoted by p and p respectively. The 
kinematic viscosity is v=p/p. These equations are known as the Navier-Stokes equa¬ 
tions (Navier 1823; Stokes 1845). They are derived on the basis of the continuum 
hypothesis (see, for example. Lamb 1975; Batchelor 1967). 

We seek the solution of (HH) as an initial-boundary value problem in the space- 
time domain denoted by ilx [0, t > 0], where Q. has a smooth impermeable boundary 
do.. The initial condition is given by 

u(x, t=0) = Uo(x) G X G n, (1.2) 

and the no-slip boundary condition is 

u(x, t)=0 X G (1.3) 

The postulation that Uq is a smooth function with compact support is to simplify 
the technicalities in the subsequent analysis. 

The vector quantity, known as the vorticity lu = (wi,a; 2 ,W 3 ), is related to the 
velocity by w = Vxu. The vorticity field must be solenoidal because of the vector 
formula V.(Vxu) = 0. The vorticity is a direct consequence of conservation of 
the angular momentum of fluid particles. Because of fluid viscosity, particles are 
subjected to shearing force during motion. Locally, they behave like a rigid body 
in rotation, at the same time, their moments of inertia redistribute as a result of 
deformation. The two effects give rise to the dynamic equations, 

— i/Awi = (w.V)ui — (u.V)a;i, f = 1,2,3. (1.4) 


Because the pressure has been eliminated, the vorticity equations are more acces¬ 
sible in deriving a priori bounds which are critical in global regularity theory of 
the Navier-Stokes equations (Lam 2013). In particular, it is verified that the total 
vorticity is an invariant in any fluid motion with smooth initial velocity of bounded 
energy. 



(wi -I- a;2 -I- W3)dx = 0. 


Consequently every vorticity component must be finite. By virtue of the Sobolev 
embedding theorem, the invariant principle can be extended to show that 


w(x,t) G C'^(x) C'g'(t), t>tc>0, 


(1.5) 
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where tc is the small time limit for smoothness in the classical regularity theory. 
Since the vorticity is a fundamental quantity in fluid motion, we need an inverse 
relation which enables us to compute the velocity and the pressure. In view of the 
incompressibility hypothesis and of the vector identity 

V X w = V(V.u) - Au, (1.6) 

an elliptic equation can be derived and used to bound Vu in space because di- 
vorticity, Vxw, is a sub-set of Jacobian, Vw. The space-regularity in the velocity 
gradients ensures the pressure gradient Vp is suitably bounded. From the momen¬ 
tum equation, we assert that 

||Vu(.,t)||^^(^3) t>0. (1.7) 

In conjunction with the classical local-in-time regularity, the global well-posedness 
follows for initial smooth data of finite energy (see, for example. Doering & Gibbon 
1995; Constantin 2001; Heywood 2007). By a similarity reduction procedure, it has 
been shown that, at high Reynolds number, the vorticity field is nothing more than 
an amalgamation of viscous shearing motions, giving rise to multitudinous spatio- 
temporal scales. In the sequel, we shall work with infinitely differentiable bounded 
functions as no finite-time singularity can develop in solution of the Navier-Stokes 
equations. Viscosity of fluids, no matter how small it may be, activates energy 
consumption; the non-linearity in the equations facilitates the dissipation by self- 
multiplying viscous shears into smaller sizes. In fact, the continuity defines a zero- 
sum of three extended reals, dui/dxi,i=l,2,3 at any given instant of time. None 
of them can go out of bound as the “measure” oo—oo=a, where a is a finite real, 
is undefined in analysis. The physics is clear; any motion due to an incompressible 
flow of finite energy cannot instigate an unbounded plenum adjacent to an infinity 
vacuum in any finite region whose flow rate of matter is finite. 

The space-time vorticity boundedness, to G L°°(x) L°°{t), holds in the limit 
of vanishing viscosity Hence the Euler equations (/r=0) cannot blow up in 

finite time according to the BKM theorem (Beale et al 1984), assuming that there 
exists a fluid of zero viscosity. The impossibility of singularities in inviscid flows of 
finite energy is consistent with Helmholtz’s vortex theorems (1858) and Kelvin’s 
circulation principle (1869). Owning to lack of energy dissipation, flow motions 
described by the Euler can only re-distribute their energy. Thus it is comprehensible 
that, at least locally, inviscid flow fields may evolve into less smooth states compared 
to their viscous counter-parts. Nevertheless, the Euler equations do not constitute 
a self-contained model for the physics of turbulence where diffusive dissipation by 
viscosity is of importance in vorticity proliferation. 

2. Velocity-vorticity formulation 

In the present paper, we deal with flows in the presence of solid surface and we 
must strictly impose the no-slip boundary condition m- We first establish that, 
in the “kinematic” system 

V X u = w, V.u = 0, u = 0 X £ do., (2-1) 
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the velocity can be recovered uniquely from 

Au = —V X (jj, (2.2) 

where the region Q, is simply-connected having a boundary dO.. The necessary 
condition to confirm (1^ is the vector identity (HSl). The continuity implies that the 
vorticity field in D. is solenoidal V.w = 0. Denote G by Green’s function satisfying 
the homogeneous boundary condition. Conversely, if we are given a bounded vector 
function A in D. The Poisson equation, Au = — VxA, is inverted 

u(x) = f G(x,x')Vx'XA(x')dx'= — / VxXG(x',x)A(x')dx', (2.3) 

Jn Ja 

where the boundary term from integration by parts vanishes because function G is 
zero on the boundary, and the last integral holds in view of the symmetric property 
of Green’s function: G(x, x') = G(x',x). Thus the incompressibility V.u = 0 is 
verified if we take divergence operation. Taking curl on the Poisson equation driven 
by A, we obtain 

A(V X u) = AA- V(V.A). 

If the vector A is solenoidal and coincides with w a.e., then the function A = uj + At, 
generally holds such that A = uj, and AAi, = 0, V.Ab = 0, and this system has no 
boundary conditions. Given Ab as a vector, we must have Vx Vx Ab = 0 or the last 
elliptic system reduces to 

V.Af, = 0, V X Ab = a, 

where a is a constant. Thus the following condition can be established: 

^f,(x £ dD.) = a dx = a vol(n), (2.4) 

Ja 

which defines a vortex sheet over the boundary. In application, we have two scenar¬ 
ios: (I) The volume £1 is unbounded. Since we are interested in fluid motions of finite 
energy, there can be no singular vortex sheets on the boundary in view of (II3| and 
(US. The function A is nothing but the vorticity. (2) For domains of finite vol(£2), 
the function Ab can be fixed by initial data (say at t=0). The problem is to trace 
the spatio-temporal flow development from a set of initial-boundary conditions by 
solving the Navier-Stokes equations. In particular, if the initial vorticity in £2 is 
zero, the system of (ED and ED is well-posed in any singly-connected region with 
smooth boundaries; the velocity can be recovered uniquely from the vorticity ED- 
Our ultimate aim is to solve the dynamic equation ED or ED in £2 where both 
u and UJ exist and evolve in real fluids. The use of the dynamic no-slip boundary 
condition is therefore justified as homogeneous Dirichlet condition (11.31) is valid at 
every given instant of time. 

In addition, imposing ED has another advantage which simply fixes the velocity 
field in the Helmholtz decomposition. Let us write 

u = Vcf -I- V X \|/. 

The function \|/ can be determined by taking curl operation so that A\|/ = —V x u = 
—UJ as long as the \|/-field is solenoidal. The continuity demands that (jj is harmonic 

A(l) = 0. 
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In it is postulated that the “potential” (j) decays at infinity, or we consider 
only finite energy solutions as a requirement for physics. Thus a harmonic function 
bounded in must be a constant by virtue of Liouville’s theorem. The contribution 
from the irrotational part drops out implicitly. For real fluids in the no-slip 
condition adds the boundary condition V(/)(x) = 0, Vx G dQ.. The solvability 
constraint for the Laplace equation, 



holds. Hence the only solution for the equation is a constant. Briefly, the formulation 
uniquely specifies the velocity or the vorticity field for real fluids /i>0. 

3. Dynamic equations in cylindrical co-ordinates 

In a cylindrical co-ordinates system (r, 0, z), we consider the initial-boundary value 
of flow evolution. The velocity and pressure (u,p)(x, t) = {u,p){r,9, z,t). Let u = 
{u, V, w) denote the velocity components in the radial, circumferential and axial 
directions respectively. Without loss of generality, we take the pipe radius (denoted 
by a) as unity. The domain of interest consists of the interior of a semi-infinite 
circular cylinder 


(3.1) 


0<r<l, O<0< 27r, 0 < z < oo. 


We denote the wall at r = 1 in 12. The cross-section at the location z = 0 is the 
pipe entry. The pipe surface is assumed to be hydraulically smooth. The choice 
of the co-ordinates attempts to model, as closely as possible, practical laboratory 
conditions. In many pipe-flow experiments, controlled water or air from a reservoir 
is forced into the pipe under external pressure differences. By experience, good 
flows may be generated in practice by maintaining a constant mass flow through 
the entire pipe test section while, at the same time, care must be taken to minimise 
fluctuations in the flow rate. The flow far downstream is inevitably in turbulent 
stage and the presence of the turbulence needs to be accounted for in theoretical 
treatment. The axial length is typically 100-1000 times of pipe diameter, depending 
on specific set-up. Although the practical pipe length may not be long enough for 
the flow to decay completely at the far end of the pipe, the flow field close to the 
inlet must be, to a good approximation, independent of downstream flow. 


Momentum equations 


The continuity equation reads 

^ 1 diru) 1 dv dw 


V.u = - 


r dr 



(3.2) 


The Navier-Stokes momentum equations are 



(3.3) 
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where A is the cylindrical Laplacian 


1_5 

dj.2 ^ Qj. ^2 QQ2 Q^2 ’ 


(3.4) 


and the differential operator A' = A—1/r^, and S/B = {drB,dgB/r,dzB) for a 
scalar B (see, for example, Batchelor 1967; Schlichting 1979; Wu et al 2006). 

The pressure p satisfies the following Poisson equation: 


Ap 

P 


nd +^[rT0 + r) 

/I du dv nN2 /du dw N2 /dv 1 9u> n 2 


(3.5) 


Hence the pressure can be determined uniquely with respect to a reference pressure 
once the velocity gradients are known. The elliptic equation contains no time-wise 
information; any variation in the velocity gradients causes instantaneously changes 
in the pressure. This apparently unphysical causality effect is a consequence of the 
continuity (13.2p . If the density of the fluid is assumed to vary with the pressure, as 
given in a state equation, the pressure variation propagates at the local speed of 
sound. In standard laboratory conditions, the speeds of sound for pure water and 
air are about 1500 m/s and 340 m/s respectively. Thus the incompressibility is a 
well-suited hypothesis in pipe flow experiments. 

We suppose that the motion starts impulsively from rest at time t = 0 and the 
flow is initiated afterward. The initial condition is taken as 


u (r, 0, z > 0, t < 0) = 0. (3.6) 

Equations (EH) and (ESI are to be solved subject to the following boundary con¬ 
ditions Vt S (0,T < 00 ]. The velocity must satisfy the no-slip condition at wall 

u (r = l,0,z > 0,t > 0) = 0, (3.7) 


and decays at infinity 


u —0 as z —>■ 00 . (3.8) 

The decay is plausible in physics for flows of finite energy. In addition, the velocity 
inside the pipe is assumed to be bounded everywhere. 


u (0 < r < 1, 0, z) I < 00 . 


(3.9) 


Because of the azimuthal symmetry, it is evident that the velocity must satisfy the 
periodic condition 


u (r, 0 = 0, z) = u (r, 0 = 27r, z), (r, z) G 12 (3.10) 

at every instant of time. For t > 0, we postulate that the flow condition at the pipe 
entry can be specified as smooth functions of position and time 

u (r,0, z = 0,t) = Ue (r,0,t) = ( We Ve We){r,6,t). (3.11) 

Since the entry velocity is also governed by the equations of motion, we require 
that V.Ue = 0 and the initial energy is finite. 
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Vorticity equations 

The vorticity components in the cylindrical co-ordinates are 

1 dw dv du dw 1 d{rv) 1 du 



r 89 dz’ ^ dz dr’^ r dr r 89 


(3.12) 


in the r, 9, z directions respectively. In the vector identity, Vxw = V(V.u)—Au, the 
quantity on the left-hand side is known as the di-vorticity which links the velocity 
and the vorticity 


(3.13) 


Au = —V xuj 


for incompressible flows. In terms of the di-vorticity, the Navier-Stokes momentum 
equations dlSl) can be re-written as 


(3.14) 


8tu -I- w X u -I- vVxuj = Vx, 


where x = ~(u^/2-|-p//9), is the Bernoulli-Euler pressure (cf. (I3.5|) l. Taking the curl 
on (13.1411 and (13.2L and making use of the vector identity V.(Vxu) = 0, we obtain 
the dynamics equations for the vorticity components 



dtC, - J^AC = {uJ.V)w - (u.V)C = Z = Z, 

where the vector components {X Y Z) = X, and {X Y Z) = X. Note that the 
continuity equation is satisfied. In addition, the pressure drops out of the vorticity 
formulation (13.1511 so that the number of unknowns is now three, instead of four. The 
dynamic pressures, 8tp, 8t8rP and the other derivatives, must be found, indirectly, 
from one of the momentum equations in (13.311 . 

As our motion starts impulsively from rest, condition (13.611 means the initial 
data for vorticity are identically zero 


oj {r,9, z > 0,t < 0) = 0. 


(3.16) 


At any subsequent instant during the evolution, the vorticity satisfies the periodic 
condition 


uj {r,9 = 0, z) = oj (r, 9 = 27r, z). 


(3.17) 


From the velocity distribution (13.1111 . the vorticity at the entry can be calculated 


V X Ue (r, 9,z = 0,t) = oje {r, 9, t), 


(3.18) 


where we have explicitly expressed Ue or Wg as a function of time. The axial gradient 
of the entry vorticity is given by 



(3.19) 


The left-hand sides in (I3.15|l are diffusion operators. During the vorticity evo¬ 


lution, the pipe wall acts as a source for vorticity in view of the no-slip condition. 
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However, the amount of the vorticity generated across the wall per unit area per 
unit time is not known and clearly depends on the dynamics of the whole flow 
field. Theoretically, the wall vorticity must constituent part of the solution. For 
convenience, we designate the vorticity boundary value as 


u} {r = l,9,z,t) = ujb{0,z,t). (3.20) 

To simplify our notations, the following short-hands are helpful: 



where the middle integral is over the entire pipe wall and it is used to specify the 
wall vorticity while the last over the inner cross-sectional area. Similarly, space-time 
integrals are 



In our presentation, we use the following notations for the independent variables 
x=(r,6»,z), y = (r,6»), z. = {0,z). 

Use of them becomes clear when the domain of integration is stated. 

In the present paper, we intend to develop an analytical theory for the flow- 
field evolution from the commencement of motion, given the prescribed data. The 
fluid dynamics in circular pipe has been subjected to intensive investigations for 
more than 130 years, starting with Reynolds’ experiments (Reynolds 1883). We 
refer the reader to a recent review and the references therein on the latest devel¬ 
opment (Mullin 2011). It is an experimental fact that pipe flow is too mysterious 
to comprehend. Little may be extracted from experiments as guides for theoretical 
treatment. It would be fair to acknowledge that laboratory experiments using vari¬ 
ous exploitation techniques are almost exhaustive as far as incompressible flows are 
concerned. It is lack of a rigorous mathematical framework which hinders progress 
because of the difficulties to reconstruct the non-linear characters of fluid dynamics 
from test data. 


4. Kinematics of vorticity induction 

The components of Poisson’s equation in (13.1311 have the explicit expressions 


A'u = - 
A'v = - 
Aw = — 


/I d( dr] 2 dv\ 

\rm “ ^ ^ M/’ 

/d^ 2 du\ 

^ ^ M/’ 

V r dr r d9 


(4.1) 


which are viewed as a system of elliptic equations subject to conditions (lT7l) on 
pipe wall and (|3.11D at pipe entry. Green’s function for the Laplacian A can be 
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found by the method of separation of variables. In view of the no-slip condition at 
the wall, it has the form of 




1 -/o(Afcr) Jo(A/e/) 




Afc Ji (Afc) 


Eo{z,z') 


oo oo 


2EE 




El k^i Crn,k (Jn+l {o’n,k)) 


cos (n(0—0')) E{z, z')j , 


(4.2) 


where 


Eo{z,z') = exp (-Afe| 2 ;-z'|) - exp (-Afc|z-b 2 ;'|), 

E{z,z') = exp (-cr„,fe|z-z'|) - exp (-cr„,fe|z-|-z'|), 

and the constant, Afe,fc = , is the fcth positive zero of Bessel function 

Jo{^) = 0, an,k the fcth positive zero of J„(cr) = 0. The zeros, cr„_fe, are larger 
than Ai = 2.404826 (see appendix A for some examples). Similarly, Green’s func¬ 
tion for operator A' is found to be 


D(x.x') = i f; g (»(«-«')) E{z.z'). (4.3) 

^ n=l fc=l <^n,k \Jn+l\0'n,k)) 

There exists a rich collection of literature on Bessel functions, see, for example, 
Watson (1944) and Olver et al (2010). Bessel functions Jn{x) are entire functions 
of the argument x, 


k=0 ^ 


(-1)'= (xV4)* 


n -I- A: -f 1) ’ 


n > 0, 


(4.4) 


and 




(-l)'=(xV4)''r(2n-b2A:-bl) 


fc=0 


(r(n-b fc-b 1))2 r(2n -b A: -b 1) ■ 


(4.5) 


The following recurrence relations are well-known (for a; > 0) 

Jn+l{x) = {2n)Jnix)/x - Jn-lix) and J^ix) = Jn-l{x) - nJn{x)/x. (4.6) 
The denominators in the Green functions are derived from the orthogonal relation 

2 / xJn(^0'n,k^^ { 'G,+ l (^n,A:) ) 5 

where 6kj is Kronecker’s symbol. For the Bessel fnnctions of all order n > 0, 


pOC 

/ Jn(x)dx = l, 

Jo 

and yn(a;) 

<1 for 0 < a; < oo 

For small arguments, we have 



Jo(0) = 1, 

Jn{x) —>■ x" as 

a; —>■ 0 n > 1. 
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The Bessel functions oscillate for large arguments 

Jn{x) ^ yj — cos (^ 2 ; — (' 2 '*” 4 )’’') a; —>■ 00 . 

By the series expansion for the Bessel function, we find 

^ I y / \\ ^ ^n,k ^ 7 \ 1 

- < —Mn-i(o-n,fca;) < ——, n > 1, fc > 1. 

X z ' '2 

We have carried out some numerical experiments for the double sum 


(4.7) 


oo oo 


EE 




For n,k ^ 0(100) and 0<r', r'<l, all the computations can be done efficiently. 
There have been no numerical difficulties. 

The inversion of the last equation in (HD) is given by 


^') (7 + ^- 7 ^) r' [d,,G] (x, y) u;e(y)dy 

= J 77 (x') - r' {dr'G) ? 7 (x') + {de'G) ^(x'))dx' + We(x), 


(4.8) 


where the functions, 

dr'G{x,x') = dGjdr', 9e'G(x, x') = dGjdO', (4.9) 

are obtained from integration by parts. The boundary terms vanish in view of the 
periodic condition for vorticity (13.171) . and the properties of the Bessel functions. 
Moreover, the notation [9z'G] = dG/dz'\^,_^, and we have introduced the notation, 

Ve(x) = (Ge K We){x), 

to account for the contribution from the velocity at the pipe entry. To continue the 
process of inversion, we reduce the first two equations in (14.11) to 

u(x) = J (^r'{d^>D) ? 7 (x') - {dg>D) C(x')^dx' -2 J u(x')dx' + Ge(x), 

u(x) = J r'[{d^>D) C(x') - (ar'I?) ^(x')^dx'+ 2 J M(x')dx'+ 14(x), 

(4.10) 

where the derivatives on Green’s function D are derived in the same manner as 
those in (HU), and use is made of the velocity periodic condition (13.101) . and the 
decay of Green’s function at z' —>■ 00 . In particular, the function dgiDjr' remains 
regular as r'—>-0 by virtue of (14. 4p and (14.7p . The induced velocities by the entry. 
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Ue and Ve, are found in an analogous manner to that by We- In (14.8L we replace 
the kernel G and We replaced by D and Ue to obtain [ 4 , 

C/e(x) = / r'[(9;,/D](x,y) Ue(y)dy, (4.11) 

dn® 

where [clz'Il] = dD /The component 14 is obtained analogously. In view 
of the incompressibility, the induced velocity Vg is completely determined and fixed 
for given geometry. It is evident that the influence of the entry data of finite energy 
is restricted to a portion of the pipe downstream of the entry and diminishes at 
large distance z —>• 00 as both Eq{z, z') — >■ 0 and E{z, z') —>■ 0. 

Because of the rotational symmetry from 0 to 27r, we readily verify that 

^^dxdx'= ^ ^ ^ de^Dddde'^dzdz'drdr'= 0. 

From the properties of the zeros for Jn{o'n,k), neither 1 nor 2 is an eigenvalue of 
the homogeneous operator A'u = 0. The equations in (14.101) are a pair of Fredholm 
integral equations of the second kind with continuous bounded kernels. They may 
be re-written in vector form for the unknown vector U = {u v) 

U + J^*U = je*uj + V = ^, (4.12) 

where V={Ue 14 )j and the ad hoc operator * stands for integrals of matrix kernel 
multiplying a vector function. (It should not be confused with the usual convolution 
operator.) The domain of integration is clear by considering the kernel and the 
function; it is over the space x. Explicitly, the kernel is a non-zero 

square matrix with zero trace, 

0 2de>D/r'\ 

\-2d0,D/r' 0 

and x') is 2x3, 

0 r'd^^n -de'D \ 

dr'D 0 r'd^>D ) ' 

Pre-multiply (I4.12p by we obtain 

Carrying out the multiplication once more, we get 

After fc-times multiplications, we write the result as 

k 

U = Y^ Jif* ■ ■■ + ^ - * ■ ■ ■ * Jf* U. (4.15) 

i fold fc+1 fold 



(4.13) 


(4.14) 
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As shown earlier, the boundedness for ds'Djr' implies ||<Ao. Since we are 
interested in fluid motions of finite energy, we consider ||ue||z,oo<i3o, then the tail 
on the right in (14.1511 vanishes as 

o 

*■■■* 1* 1 ) < —2— —^0 as fc —>■ oo, 

fe+l fold 

where I is the identity matrix. By the Sobolev embedding theorem, we may make 
reference to the L^-theory for Fredholm integral equations of the second kind (see, 
for example. Chapter 2 of Tricomi 1957). The operation sum must converge in the 
limit of /c —>■ oo. Let us denote the infinite sum of the multiple operations by 
and it is called the resolvent kernel of Jif. For bounded vorticity, the solution of 
the integral equations (14.121) or (14.101) is given by 

U = ,je*uj + Jf)*uj + V + W* V. 

Combining the result with (14.811 , we write the relation between the velocity and the 
vorticity asu = 'F*uj + 'W* Vg, or 



This relation holds for every instant of time but it does not contain time informa¬ 
tion since equation dUD is elliptic. This is a consequence of the incompressibility 
hypothesis (1321). We draw our attention to the fact that the vortices, ^ and r], do 
contribute to their velocities u and v in contrast to the velocity induction in the 
Cartesian co-ordinates. The reason is that u Laplacian is linked to the azimuthal 
gradient in v, which is partly driven by vorticity component When a unit-mass 
fluid particle at r = ro from the pipe centre is rotating in the r—6 plane instan¬ 
taneously, the term jr^ in the u-momentum equation in (13.1411 is the centrifugal 
acceleration on the particle toward the centre. 

It may have been a generalisation of the potential theory or the electromag¬ 
netism in classical physics, suggestion has been put forward to represent the veloc¬ 
ity field by a solenoidal vector stream function, where u = V x T. Instead of (14.11) . 
we have 

A'F(x) = -V X u = -w. (4.17) 

Although this set of lower-order equations appears to be simpler, its general use 
can be trick and hence is not recommended except for flows in unbounded space . 
The reason is that no quantitative knowledge of the stream function on solid walls 
is available; the inversion of (14.1711 cannot be readily obtained. The integral of u 
contains arbitrary constants, depending on local flow condition. On the other hand, 
imposing the no-slip condition involves two Green’s functions for every component 
of ‘F and hence the evaluation of these Green functions becomes substantially more 
complicated. 
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5. Vorticity due to entry flow 

Because of the specific characters of our co-ordinates, we prefer to write Green’s 
function for the diffusion kernels in (13.151) as 


H(x,x',t) = Zi{z,z',t) r' K{r,r',6,9',t), 
where H is diagonal, and we denote its components in a row matrix 

n = {H H H 3 ). 

The function, 

1 r / ( {z + z'f 


Zi{z,z',t) = 


y/Hdir 


I exp 


4i/t 


) + exp 


Avt 


(5.1) 


)}, (5.2) 


is the gradient of the fundamental solution of heat equation in one space dimension 
{z > 0). The second part of the Green function satisfying the Neumann homoge¬ 
neous condition on the wall is written in the diagonal matrix 


K=iK K Ks). 

The first two elements are K = I/tt Kq, where 


00 00 






cos (n(6»-6»')) exp i/f), (5.3) 


and 


11 .^ Jo{l3mr)Jo{l3kr') , ^ , 

^3 = - + - > \ -^ exp (-/3^ vt) + 2 Ko. 

TT TT ( ^ I ^ i H \ 


m—1 


JliPm) 


(5.4) 


The constant, (3m, m = 1, 2, • • • , is the mth positive zero of Bessel function J'q{(3) = 
0, and an,m the mth positive zero of Jn{a) = 0. The zeros, an^m, are larger than 
Pi = 1.841184 (see appendix A). 

As a first approximation, we investigate the vorticity field caused by the entry 
vorticity (13.161) and, for the time being, we ignore the contributions from the wall. 
Then the dynamics equations (13.151) can be expressed in integral form: 


w(°)(x,t) = 


[ [ Tie{x,y,t-t')ve{y,t')dydt' 
Jo Jd.<Si 


(5.5) 


+ H(x, x',t—t')X^°^(x',t')dx'dt', 
Jo Jn 


where Hg is a diagonal matrix. We denote the elements by 

He = (A A A3), 


and 


ne = -ur'K[Zi] r' K(r, /, 9,9', t - t') Zo(z, t-t'), 

where t > t' ^ 
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As implied in the last integral term represents the non-linear 

interaction in the vorticity field. 

By contrast with the vorticity dynamics in the Cartesian co-ordinates, two extra 
terms arise in (I3.15L namely, 

We expect that they are invariant with respect to the rotational symmetry. For 
instance, we obtain the following result by integration by parts: 



/a® 


/ATfr, r',6»,6»')4(r',6»')dr'd6»' 


(5.6) 


for fixed time t'. We draw our attention to the fact that this reduction is performed 
independent of the heat kernel Z\. The first two vorticity components in can 
be reduced to integral equations 


f) = 2z/ / f x', <—t') t')dx'dt'-I-A */e-I- 

Jo Jn ’’’ 


77 ^°^(x, t) = — 2z/ / [ Zi ^ (x, x', t—t') t')dx'dt' + A* + JI* 

Jo Jn r' 

(5.7) 

This is a system of two Volterra-Fredholm integral equations of the second kind (cf. 
(I4.12[) l for the unknowns over space-time. We note that Kjr' is regular 

at r' = 0 because Ji oc r near r = 0. The operator * indicates the integration over 
space-time (x, t). We express the kernel of the system as 

0 2vZide,Klr' 

-2vZi dg^K/r' 0 

Evidently, the kernel is a continuous bounded function of x and t > 0. Let K be 
the resolvent kernel of fco- On the basis of standard theory for parabolic differential 
equations (see, for example. Chapters 1 & 2 of Friedman 1964), we find that 

\iJ\ ^ ( iz±z')^\ 

I I (i/t - lyf )i/2 I “ - t')) 

for some constant C, and 0 < k < 1. Denote 

Ri = K * A, R2=K * H. (5.8) 

Then we verify that 

ISI, C, / 

I'*'! < MiV?)j- 
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and 

1^2 I < C2{l^t - Vt' exp - K2 

where Ci and C 2 are constant, and 0 < ki,K 2 < 1- The zeros of the Bessel func¬ 
tions, an,m, are large positive values. Thus the term in Green’s functions K has 
exponential damping effects in time. The pair of integral equations can be solved 
now. The result is combined with the ^-component to give 

^(0) = A * /e -h .Rl * /e + ^1 * ffe + + ^2 * + ^2 * , 

=24*ge + ^l*3e + ^l*/e + ^f*T(°) (5-9) 

= A3*he+i?3*^^°\ 

where the entry vorticity in the radial and the azimuthal directions is filtered by 
the Volterra-Fredholm operator. For example, the second term in the first equation 
on the right is, for t > t' > t", 


i?i */e(x,t) = / / K{x,x',t-t')A{x',y,t'-t'')fe{y,t")dydx.'dt'dt" 

Jq.t Jnf 


= / ^ ^i(x,x',t-t')/e(y,t')dydt', 

Jnf 

where we have interchanged the orders of integrals, assuming Fubini’s theorem. 
Similarly, the fifth term is 

R2*X^^\x,t)= ( f K{x^x',t—t')H{x\x" a', t")dix" dt" dx'dt' 

J Q.'j- J Q,j- 

= f R 2 {x,x', t—t')X^°\x',t')dx'dt'. 

J Q.'j^ 

We write the system (15.91) in the short-hand form 

w*- 

where the matrix kernels are 


G, = 


and 


G, = 


— Ge * Ve 

+ Gg * 

xw, 

(5.10) 

A d- R\ 

Ri^ 



Ri 

A -j- Ri 

0 

, (5.11) 

0 

0 

kl3 / 


H + R 2 

i?2 

0 ^ 


i?2 

H + R 2 

0 

(5.12) 

0 

0 

) 

I 


From result (15.6p . definition (EH) and properties of kernels A and H, we see that 
the elements of Gg and G^ are proportional to 

exp (- and exp (- (5.13) 
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respectively. To save space and typing, we write 

as (u, w)^^\ 

Let us consider component in (15.91) . First 


H 


t) = f iL(x,x',t—t')^(a;.V)u — (u.V)^^ (x', t')dx'dt' 

= [ Vx'i?(C'c + —Cw)^°^dx'dt', 

Jq.t 


where we have performed integrating by parts in the space variable, and made use 
of the no-slip boundary condition for u, the periodic condition for u and w, zero 
boundary value for K at r'=0, and the solenoidal conditions V.u=V.u;=0. Note that 
the two terms from the vorticity convection and stretching cancel each other. 
By the same token, we find similar expressions for the other two components: 

iJ * F(°)(x, t) = f -b ryu — (u — ^u)^°^dx'dt', 

J 

Hs* = ( Vx'iL3(CM-b C'l'— — ?yw)^°^dx'dt'. 

J Q.'j^ 

Second, consider 


+r(°))(x,t) = [ i? 2 (x,x',t-t')((w.V)u-(u.V)e 
Jqt ^ 

\ ( 0 ) 

+ (a;.V)u — (u.V)ry + {f]u — ^v)/r'j (x', t')dx'dt'. 

Carrying out similar expansion, the right-hand side is simplified as 

f +'n)w - - Cv) dx'dt'-b f ^(vu-^v) 

J J ^ 


dx df . 


There have been some cancellations among the terms arising from the vorticity 
convection and stretching. We notice that the same reduction can be made for 77 *^°^ 
component. Lastly, component is simplified as 

( 0 ) 


J Vx'Lfs (Cu + C'^ ~ (? + dx'dt'. 


Thus the system (15.101) can be separated into the individual contributions from 
the vorticity and the velocity: 

a;(0) = Ge*Ue + Ci*(F^°) ) -b C 2 * (F^°^ u(°)) -b C 3 * (F^°^ u(°)), (5.14) 

where 


F^°) = 


-(^ + 0 c e 

V -(C + C) V 

c c -iC + v) 


(0) 


(5.15) 
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p(0) 


-C -C iC + v) 


(0) 


■2 = I -C -C + 


0 0 


0 



(5.16) 


In equation (I5.14L the brackets are used to emphasise our conventions that the 
vorticity and velocity matrix multiplications are performed first. The kernels C’s 
are diagonal, and we summarise their elements in 


Ci = (Vx'i7 Vx'i7 Vx'i73), 
C2 = (Vx'i?2 Vx'i?2 0), 

C3 = (i?2/r R2/r 0). 


(5.17) 


In view of bounds (15.131) and the properties of the heat kernel Z\ in Green’s function, 
the time-wise dependence of these elements becomes clear. In particular, we have 


Vx'i? oc {vt — vt ') 


- Id’ 1 - 3/2 


exp 


{z±z'f 
4iy{t - t') 


)■ 


(5.18) 


By virtue of the inverse relation (j4.I6ll . we substitute vorticity for 
After straightforward algebraic manipulations, we find that system (|5.14ll can be 
expressed as 

u!^^l=Ge*Ve + J ^Rq (ug)-f Ri (w^®^ )^ (x, x', f—f') (x', t')dx'dt', (5.19) 

where the elements of Ro and Ri can be found from those C’s, Vij and Wij. By 
the estimate (I5.I3|1 . the term associated with the entry vorticity is found to behave 
like diffusion in one space dimension, 


Ge * 1 


,{x,t) (x{iyt) ^/2exp(^-^y 


The flow evolution really occurs in three space dimensions because r and 0 vorticity 
components, ^ and t], are closely coupled. In terms of the resolvent kernel of Rq, 
this integral equation is further reduced to 

= Ge ★ Ue-b f Q(x, x', t—t'; 0 ;^°^(x', f')dx'dt'. (5.20) 
J 0 ,- 2 ^ 

The last integral is quadratic in vorticity components, The kernels Ge and Q 
are non-linear functions of the entry velocity Ue. 

Following our treatment of vorticity development in K3 (see Chapters 7 and 8 
of Lam 2013), we notice that every integral component in ()5.20l) may be expressed 
in the form of 

3 

f f ^Aa{x,x\y,t-t')ujl°\y,t')uj^°'>{x',t')dydx.'dt', f = l,2,3. 

J J ^ 

In view of (I5.18|) . we have 

2 

Ge * Ue(x, t) ( = ■n7(x, t)) OC {vt)~^/'^ exp . 
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Similarly, coefficient matrix Aij has this identical temporal dependence. Thus the 
non-linear integral equation (15.2011 has the component form 

r r ^ (5.21) 

+ / Aii{x, x', x", (x", (x', t')dx"dx'dt'. 

JcitJo. 

These expressioirs imply an algorithm for similarity reductions; over a small time 
interval in At € [0, t], we multiply the left-hand by An and for every z = 1, 2, 3. 

The sum of the results is equivalent to the difference of and the initial data 

. The reductions can be performed as many times as we wish. Written in vector 
form, the non-linearity in the integral equation is effectively transformed into an 
infinite sum: 

t) = •n7(x, t)-t / B^°^(x, x',t—t')a;^°^(x',t')dx'dt'-I-(/(x, t), (5.22) 

J 

where q is an infinite sequence 


OO 

m— 2 

and stands for m-fold integral of An and zu. Let be the resolvent kernel 
for The unknown vorticity in equation (I5.22p can be recast and expressed as 

a;^°^(x, t) = 7 (x, t)-I-g(x, t)-I- f {x,x'+ q){x',t')dx'dt', (5.23) 

J Q,ip 

where 7 is the Volterra-Fredholm filtered entry data, and is given by 

7 (x, t) = ■n 7 (x, t)-f f B^°^(x, x',t—t')TO(x',t')dx'dt'. (5-24) 

J 


By the method of successive approximation, equation (15.231) can be completely 
solved for any given initial data at pipe entry Ve- We express the solution in terms 
of the vorticity integral powers fo[ 7 ]* = fo[ 7 (x, t)]^ in the following convergent 
series: 


(x, t) = 7 (x, t) -b 2 7 ]2 -H 10 fo[ 7 -b 62 7 -k 430 fo[ 7 

-b 3194 fo[7 ]®-b 24850 fo[7]V .. 


(5.25) 


It is convenient to write the solution as 




(5.26) 


6. Flow evolution over space-time 

After a small time t = tg, the flow field defined by (1531) or given by local solution 
(15.261) occupies the entire interior of the pipe. For convenience, we may re-define the 
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start at time t = tg. Applying Duhamel’s principle, the dynamic equations (Id. 1511 
can be expressed in integral form (t > t'): 


a;(x,t)= f f lle{-K,y,t-t')veiy,t')dydt' + / H(x,x',t)a;s(x')dx' 

Jo Jn<s> Jn 

+ f f VLb{x,z,t-t')(pb{z,t')dzdt' 

Jo JdQ. 

+ H(x, x', t-<')X(x', t')dx'dt' ( = le + /. + A + $) , t > t'. 


( 6 . 1 ) 


This equation comprises all the vorticity contributions of the initial-boundary value 
problem (cf. (15.51) 1. 


(a) Initial vorticity 

As in many physical problems where parabolic differential equations are appli¬ 
cable, the second term on the right is the contribution from the initial vorticity. 
If the interior fluid is set into motion impulsively from rest, instead of the void 
vorticity (13.161) . we must specify the initial distribution over the whole pipe wall in 
the form of a vortex sheet. The strength of the sheet is difficult to quantify without 
approximation. By considering the local-in-time shears on the wall due to the entry 
vorticity (15.251) . we prescribe the initial data for vorticity as 

a;s(r,6l, z > 0) = ( Vs Cs )ir,0, z > 0) = {r,0, z > 0,t = tg). (6.2) 

An alternative, which is popular in numerical computations and circumvents the 
uncertainty in the initial data, is to perturb an existing Hagen-Poiseuille parabolic 
profile over a part of the pipe length. Written in dimensionless form, parabolic 
distribution, w{r)=l—r'^, is widely used. The mathematical model is periodically 
extended over the whole length of the pipe. The two underlying assumptions in 
this model are that Hagen-Poiseuille flow is absolutely steady in the complete pipe 
and the parabolic profile exists at arbitrary Reynolds number. However, as an 
initial-boundary value problem, a fluid motion may start from an arbitrary initial 
profile. The dependence of the initial condition cannot be ignored as the initial 
data dominates the subsequent flow development. In fact, both the assumptions 
can hardly be justified (see m- 


(&) Wall vorticity 

The third integral states the contribution of the vorticity production from the 
boundary; it is calculated by means of time-dependent integral relation 



H{,(x, z, t—t')ipb{z, t')dzdt', 


(6.3) 


where the vorticity derivative in the r-direction at wall is evaluated from 


I>b{z,t) 


doj{x,t) 
dr r=i 


(6.4) 
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and the kernel is given by 


Hfa = i/Zi [r'K]r'=i = vZi{z,z',t) K'(r, 0, 2 , 0', x', <). 


The diagonal kernel H;, has elements 


(6.5) 


Hb = (B B Bs). 

Similarly, function K' = ( if' K' if 3 ), where if' = I/tt + Kq, and 


^ CX3 CXD 2 T / ^ 

1 X ^ X ^ ^n,m ^n\^n,m‘^) 




n=lm=l Jn{an,m) 


COS (n(0-0')) exp vt), 


and 


^3 


1 _ Jo{Pmr) 


exp {-Pi, vt) + 2if(, 


The analytic form of the kernel Green function indicates that the wall layer is 
determined by every part of the flow though the influence at the far end of the pipe 
is clearly insignificant, specifically for small data. For initial flow of arbitrary size, 
our method of solution takes the global characters into account since the vorticity 
field so solved {Lpb) contains all the interaction among the shears of different scales 
over the complete pipe length. In particular, the wall term retains the essential 
structure of the vorticity solution (15.251) . The dependence of the kernel H;, on 
viscosity v shows that it is absolutely essential to consider the viscous effect. On 
the other hand, the wall vorticity gradient illustrates the fact that the pipe flow 
field is dynamic in nature. Modelling pipe flow in terms of steady Hagen-Poiseuille 
flow of parabolic profile ought to be a crude approximation for the study of the 
dynamics; it would be mis-leading to consider the axial symmetrical flow over a 
significant part of the pipe length, particularly over Reynolds number range of 
turbulence transition. 

We refer the flow structure at the solid surface as wall layer in contrast to 
boundary layer as the latter suggests that Prandtl’s boundary layer approximations 
are involved. Evidently, the flow is everywhere fully viscous and there does not 
exist a clear-cut “free-stream” where the flow may be considered as uniform and 
inviscid. Equation dSS]) implies that the pressure must be strongly dependent on 
the instantaneous values of the velocity gradients. The shears themselves evolve in 
the transport process of convection and diffusion. 

One particular feature of the wall vorticity is rather interesting. Suppose that 
the motion is set up impulsively and we then maintain a sustained supply of the 
constant flow rate through the entry. The local vorticity diffuses out of the wall 
into the interior in the same way as heat is released from a surface. The velocity 
gradient normal to the wall contains a natural decay, i.e., the factor exp(—^ vt). 
For fluids of small viscosity, the decay near the entry appears to be time-insensitive 
and the diffusion may proceed linearly even in the presence of the wall. This is of 
course a mis-interpretation because the main vorticity dynamics is determined by 
Lpb which is non-linear and unsteady. In particular, if vorticity a;(x, t) has a complex 
structure like turbulence, the flow in the wall layer can be readily understood even 
though experimental measurements may be difficult to conduct. 
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Solution by time-marching 

We solve the vorticity evolution (EH) by considering the following iterations: 
For j = 1,2,3, •••, 


= Ie{x,t) + /p + 4(x,t) + 

where the initial data are given by 

/p“^^(x, t) = f H(x, x',t)a;^'^“^^(x')dx'. 

Jn 

The wall vorticity gradient EH is written as 


/f,(x,t)= f f Hf,(x,z,t-t')(/Jb(z,t')dzdt', 
Jo Jdn 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


where the starting value may be taken from the result of 

{9,z,t) = {pb Qb rb){0,z,t). 


( 0 ) duj^°'> 


r—1 


Solution 

The new vorticity evolves according to 

a;^^^(x, t) =/e(x, t) + 4°^(x, t) + / Ht,(x, z,t-t')(/3f,(z,t')dzdt' 

JdQT 

+ f H(x, x',t—t')X^^^(x',<')dx'dt'. 
J Q.j’ 

The updated pair of integral equations EH now read 

^^^\x,t)=2iy f f Zi ^ (x, X.', t—t') (x', t')dx'dt' 

Jo Jn 

A-k fg H -k -\- B k Pb -\- H k X^J\ 

de'K 


(6.9) 


( 6 . 10 ) 


ri^^\x,t) = —‘Jiy I I Zi ^ ^ (x, x^,t—t') ^^^^(x',t')dx'dt' 

+ Akg^ + Hkr]i°^ + B k Qb + H kY‘^^\ 


0 Jn 


Evidently, the reduction of the integral involving can be carried out in the 
same manner as we did for X*^°) (cf. (15.51) 1. There is no point to repeat the detailed 
procedures. Equation (15.141) takes the revised form of 


do) 


— Ge k Ve J- Gg * + Gft * 

+ Cl * (Ff ^ u(i)) +C2k (F^ u(i)) + C 3 * (F^ u(i)), 
where Gg is given by (I5.12L and 


( 6 . 11 ) 


B + Kb Kb^ 0 

Gh = I Kb B + Kb 0 

0 0 B3 


( 6 . 12 ) 
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and Kb = K -k B. The matrices, Fi, F2 and F3, are still given by (15.1511 to (I5.16|) 
with the current value of vorticity Since 

2 

B, B 3 oc exp 

we calculate that 

K, « (rf - „t' )■/“ exp ( - 
for 0 < K < 1. Equation (|5.19p is revised as 


(x, t) =Ge * ue + Gs * 4 °^ + Gb * 


Fi 


( 0 ) 


+ J ^Ro(ue) + (x,x',t—t') a;^^^(x',t')dx'dt'. 


(6.13) 


Although we have not yet had the full structure of the wall layer, we may infer from 
this result that it is dominated by the global vorticity field. Essentially, the last two 
integrals contain the quadratic non-linearity. The analogous equation to (15.2011 has 
the expression 


UJ 


W = Ge * Ue + G, * -t- Gb * ¥>1' 


(0) 


-b f Q(x,x',t—t';a;^^^)a;^^^(x',t')dx'dt'. 

J Q.'j' 


(6.14) 


In view of the velocity-vorticity relation (14.1611 . the non-linear integral equation for 
the vorticity can be derived: 


4 \x,t) = 


+ 


j f x', x", t')dx"dx'dt 

t/</^_ 2 


(6.15) 


where it is clear that the revised zu, 

= zu{x,t) + Gs + Gb* 

takes into account the initial vorticity and the vorticity production at wall. Accord¬ 
ingly, the function 7 can be updated to 

+ f (6.16) 

J Q,rj^ 

The solution of (16.1511 is given by 

a;W(x,t)=J^(7W), 

where the subscript refers to . Next we perturb the wall vorticity gradient 
by a small amount to 

(1) (0) I X 

ifi' =ipi' + d(pb. 
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From the dependence of the wall vorticity given in we have 
~ f ~ ~ 

—^(x, <) Ri Gb(x, t) + / B^^^(x, x',t—t')Gt,(x',t')dx'dt'< 00. 

JaT 


Thus the starting wall vorticity 


( 0 ) 


is replaced by 



dr 


{e,z,t) 

r—1 


{Pb Qb rb){0,z,t). 


With the updated wall vorticity gradient, we solve equation (16.91) to obtain the 
updated solution This is nothing more than a Newton-Raphson procedure. 
We continue our iterations until 


max I I < 

(x,t) 


fc = 1,2,3,- 


where e is a prescribed tolerance at fc = M. The final vorticity solution is denoted 
by = w(b(x,t). 


Solution 

Once has been determined, we march into time t St. By shifting the time 
origin, we carry out the similar calculations for First we update in (16.91) 
to . In view of the Newton-Raphson iteration, the wall vorticity and hence 
can be found. The solution may be written as uj^^\x,t) = .^{ 7^^^ ). To continue 
the time-marching, suppose that we have the solution for 

(^(1-1) (x,i) = jSf( yb-l) 

The evolutional solution, is determined from the current equation (16.61) : 


,U) 


(x,t) = Je(x,t)-I-Jp b(x,t)-|- f H{,(x,Z,t-t')(/?P ^\z,t')d2.dt' 

JdQ.T 


-I- f H(x,x', t—t')X*''^px', t')dx'dt', 

J 


(6.17) 


( 7 — 1 ) 

where If 
is given by 


is the vorticity initialisation from wb and the starting wall value 


Q-_i) duj^^ b 

n = 


dr 


r—1 


The first two functions on the right in 16.171) are known. The third one may be 
determined by iterations. Now equation (16.lip becomes 


=Ge*Ve + G, *a;p-b + Gb* 

-f Cl * (Fp^ +C2* (Fp^ ubp -t C3 * (Fp^ ubp. 


(6.18) 
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Since the structure of this non-linear equation is identical to (16.1111 for the case 
j=l, in view of the velocity-vorticity link (I4.16L we convert it into the vorticity 
quadratic form, i.e., 

/* /* _ ^, (6.19) 

+ j / A^p (x, x', x", (x", (x', t')dx"dx'dl' 

JnTJn 

(cf. (15.211) 1. Thus the function can be updated: 

7 ^'^^(x, t) = ■cc7^'^“^^(x, t) -I- f B^'^^(x, x',t,t')tn^'^“^^(x',t')dx'dt'. (6.20) 

«/ 

We iterate on the wall vorticity gradient using the Newton-Raphson method to 
obtain By the method of similarity reduction and successive approximations, 
we assert that the solution of (16.6p is given by 

), j = l,2,3,---. 

The solution is unique and regular for any t>0. 

7. Gradated flow regimes 

Mathematical description of turbulence 

We have demonstrated that the vorticity evolutional dynamics governed by the 
non-linear integral equation (I5TD is globally well-posed. Its solution can be repre¬ 
sented in the series 

w(x,t)= 7 (x,t) -I- 2 1 /( 7 ]^ -I- 10 V'[ 7 ]^ -(-62 V'[ 7 ]"‘ -I- 430 R[ 7 ]® 

-b 3194 F[ 7 -b 24850 R[ 7 -b 199910 R[ 7 ]® (7.1) 

-b 1649350 R[ 7 ]® -b 13879538 7 ]^° -b ., 

where 7 = ^(ue). Function Ve or We is the vorticity value at pipe entry and may be 
called the input vorticity as it completely specifies the flow evolution. 

A vorticity power term, V[ 7 ]^=V'[ 7 ]^(x,t), generally represents the fc-fold 
integral convolution on the initial data 7 . Over space-time, the convolution effect 
cultivates a “vortice”; the higher k, the smaller the vorticity scale because of the 
heat kernel Zi. The coefficient integers form an integer sequence and the integral 
convolutions obey certain combinatorics rule^. For every additional term in the 
series, the number of vortices is increased roughly by an order of magnitude com¬ 
pared to the preceding scale. The solution dm designates turbulence in the sense 
that the vorticity field is composed of interacting shears in numerous scales. The 
solution can be computed to arbitrary degree of accuracy, at least in principle. 

t This is sequence A107841 in The On-Line Encyclopaedia of Integer Sequences www.oeis.org 
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The numerical values obtained at specific space-time locations represent the mean 
quantities in fluid motion. In reality, physics demands a demarcation on the mini¬ 
mum allowable scales. It is a fact from experiment that turbulent motions observed 
in nature and laboratory can be described by ensemble-averages and statistically 
invariant fluctuations, 

u(x, t) = u(x, t) -I- u'(x, t), u'= 0, 


where an overbear denotes the mean or the ensemble-average. Both properties are 
repeatable in practice as long as the experimental samples are sufficiently large so 
that the theory of large numbers is applicable. The missing random characters in 
real turbulence is out of the question in the current continuum approach. Because 
of the viscous diffusion, we may examine the origin of the randomness in terms of 
the kinetic theory of gases. The law of energy conservation dictates the loss of the 
kinetic energy via viscous dissipation equals the increase in the internal energy of 
fluid’s molecules. The diffusive process on the shears of dissipative scale must be 
irregular because the thermal energy of the molecules fluctuates randomly (Lam 
2013). From application point of view, we notice that the structure of the wall 
vorticity ipb resembles many features in turbulent flow. 

To preserve the generality of solution (EU, we shall not, at this stage of analysis, 
classify flows into particular categories, such as rolls, streaks, hairpins and station¬ 
ary or travelling waves though these dynamic flow features are basic building blocks 
for complex flow fields at large Reynolds number. The initiation, sustainment and 
diversification of these definite details must be history-dependent and scale-specific 
at least during the early phase of flow evolution. 

Given entry data (13.111) , the magnitude and distribution of the input vorticity uie 
determines the subsequent dynamic flow development. In general, a characteristic 
velocity at every axial location throughout the pipe prior to transition must be the 
function of entry velocity Ug. As the local flow is in evolution, the local velocity is 
evolving and hence cannot be constant. In the present paper, we define the Reynolds 
number as 



(7.2) 


V 

where w = ?Zj(ue) is the mean velocity, and d the pipe diameter. This is a textbook 
definition for Hagen-Poiseuille parabolic profile. 


W/Wmax = 1 - r^, 


(7.3) 


and w = ?nniax/2. In practice, it is extremely difficult, if ever possible, to generate 
and maintain the parabolic velocity from pipe entry to end. Most pipe flow exper¬ 
iments are carried out in the condition of constant stream-wise pressure gradient 
or of constant mass-flow, for flow starting from specific inlet design to minimise 
the local disturbances. Velocity measurements in the vicinity of pipe entry indicate 
that the local profiles are not self-similar. A Hagen-Poiseuille parabolic profile or 
the fully developed pipe flow may only be realised at downstream distance of order 
of 100—1000 diameters for small initial data or lower to moderate Re. 

We would like to emphasise the fact that the similarity rule implied in relation 
(17.211 is valid strictly for the given distribution Ug. Care is needed in the interpreta¬ 
tion of conclusions drawn from different experimental measurements, particularly 
for those on the nature of laminar-turbulent transition. 
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(a) Uni-directional approximation for entry flow 

The velocity components in the radial and circumferential directions are as¬ 
sumed as identically zero over a small time interval t £ [0,Te]. The axial compo¬ 
nent is treated as a function of the variables r and t only, i.e., w = w{r,t), and 
the pressure varies with the axial variable, p = p{z,t). The differential equation is 
much simplified as 

w(l) = 0, w{r,0) = We{r), (7.4) 

where Pq = —{dpldz)lp, denotes the axial pressure gradient. The solution of the 
initial-boundary value problem can be expressed in terms of the Fourier-Bessel 
series (cf. Lam 2015), 


{r,t) = [ We{r')Go{r,r',t)dr'V f f Po{z,T)Go{r,r',t-T)dr'dT, (7.5) 
Jo Jo Jo 


where Green’s function is given by 


Go(r, r', t) = 2r' ^ exp ( - \\vt) 

n—1 


<70 (A^ir) To (Ay^r ) 




(7.6) 


This is the one-dimension version of (I5.3I) - (I5.4I) . Consider the case of constant pres¬ 
sure gradient and the velocity is normalised by Po/4. The scaled velocity is com¬ 
puted from 


w*{r,t)= [ wl{r')Go{r,r',t)dr'd-(ii't) f ( Go(r, r', t(l—T))dr'dr 
Jo Jo Jo 


(7.7) 


= Wi WG- 


(If the initial data We=0 or the motion is assumed to start impulsively from rest, see 
Szymahski 1932 and §4.3 of Batchelor 1967, the second term produces a parabolic 
profile in the limit of vt^oo, i.e., mmax=l in the distribution (ESI).) In figures [Hand 
m we present the time evolution of two constant entry profiles. Briefly, the actual 
long-time evolution depends on the relative strength of the two contributions; the 
first term involving initial data w* may well be dominant in the solution (El) and 
hence renders the parabolic profile irrelevant. 

To relate our theory of uni-directional flow ED in the vicinity of pipe entrance 
to experiments, we set 

nt ^expt 1 ^eff 1 

a Rr d Re' 

where the quantity ccefi is an effective distance measured from an origin of zero 
viscous thickness at the theoretical “entrance”, and Rr is the Reynolds number 
based on pipe radius. It is known from experimental work (Smith I960; Pfenninger 
1961; Tritton 1976; Mullin 2011) that pipe inlets are streamlined designs to minimise 
irregularity in starting flow. The present approximate theory does not model any 
specific designs and hence the velocity profile is assumed to start some distance 
upstream of the entrance. At entry location, z=0, the local velocity then satisfies 
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w = 0.5 



w, 




w 


G 


w 


Figure 1. Influence of initial data on flow regime near pipe entry by theory of uni-directional 
flow. For scaled entry velocity: constant distributions = 0.5 and = 4.0. The pipe 
velocity is given by (TTTl) for several fixed values nt. As time nt increases, the initial velocity 
decays and the wall layer grows in thickness while the flow driven by the constant pressure 
gradient dp/dz builds up. The combined profile w* is shown in the last row. As vt^oo, 
the initial influence diminishes and the pressure gradient becomes dominant until the total 
profile approaches to the Hagen-Poiseuille parabolic distribution (ESJ. For small initial 
data, our solution explains why the fully developed pipe flow at large ut has a parabolic 
velocity distribution. The developed flow may be treated as history-independent after 
sufficiently long time. 


the no-slip condition. This practice results in a shift in quantity vtjo?. Comparisons 
in figures [3] and 0] show that the evolutional solution (17.711 is very satisfactory over 
a wide range of vt despite of the use of flow conditioners in the experiments. The 
initial-boundary value formulation of the Navier-Stokes equations must be well- 
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w = 4.0 



Figure 2. Importance of the initial unsteady phase in pipe flow. The magnitude of entry 
velocity is increased to wt = 4.0 so that the developing pipe profiles are far from parabolic 
over almost all time. For flow development at high Reynolds number, whether the local 
flow is laminar or turbulent depends strongly on the initial data. The Hagen-Poiseuille 
proHle (ES) is just a particular case of developing flows, that can be generated in long 
circular pipes. 


founded. Even in the present approximation, the solution is an excellent description 
of the initial pipe flow development over substantial entry-length. 

Our chief interest here is to study how turbulence transition initiates naturally. 
We wish to elucidate the way in which a streamlined flow structure evolves into a 
complex vorticity field in the absence of any disturbance. The experiments by Durst 
& Uensal (2006) are a synopsis which demonstrates many important aspects of the 
natural transition in well-crafted test ambience (figure [S]). 
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(a) R =0.01 (b) R =0.02 



(c) R =0.03 (d) R =0.04 



Figure 3. Theory of uni-directional flow for initial evolution in circular pipes. Experi¬ 
mental results are averaged values of the data given in Smith (1960) at Rr = wajv 
= 0.01, 0.02, 0.03, 0.04. Note that the data were taken from three collections of experi¬ 
ments, where the test conditions were very different, and, in particular, the detailed inlet 
geometries were not fully reported while the initiation flows were specifically tailored. In 
the theory, the entry velocity is considered as a constant value wj = 0.6. The full lines 
are calculations at vt/a^ = 0.015, 0.025, 0.035, 0.045 for a pipe of radius unity with 
sharp-edged inlet. The broken blue line is solution at vt/a^ — 0.25 which indicates a pos¬ 
sible decay of the inlet flow. As i/t —>■ oo, the flow can be treated as fully developed and 
tends to parabolic profile (17.311 (broken black line). The two broken lines are given for 
reference only as the uni-directional flow is supposed to exist in a small time interval from 
the start of motion. 


Distinct flow stages 

In his pipe experiments (Reynolds 1883), Reynolds established that there are 
two fundamental categories of flows: laminar and turbulent. Laminar flows are char¬ 
acterised by well-defined streamlines. Turbulence is observed on macroscopic scale 
as a fluctuating random motion consisting of various flow scales. The intermediate 
flow stage between these two flow states during flow evolution is known as transi¬ 
tion. The process of transition from the laminar state to turbulence attracts most 
attention. The reverse scenario is known as re-laminarisation. On intuition ground, 
turbulence transition must fusion the streamlined structure and the growing gra¬ 
dations of finer-scale flows. As established from experiments, transition appears to 
initiate in localised regions in intermittent manners. There are patches of irregu¬ 
larly fluctuating flows during evolution; these are commonly referred as turbulent 
puffs or slugs or spots which are often recognised in wall-bounded shear flows. 

For a flow of small Reynolds number or a fluid motion with small initial data, 
the vorticity distribution is given by the leading term 7 in (IZ3D, and the flow is 
said to be laminar. The dissipation of the kinetic energy is dominated by diffusion 
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(c) R =0.03 (d) R =0.04 



Figure 4. Theory of uni-directional flow for initial evolution in circular pipes. See previous 
plot for experimental data and theoretical computations. The entry velocity is a boundary 
layer w* = 0.6(1 — Compared to the case of constant velocity, the differences are 

negligible, as both theoretical profiles are assumed to start somewhere upstream of pipe 
entry, the starting profiles at z=0 ought to be almost identical. From the experiment point 
of view, it is much easier to produce a constant velocity in the vicinity of pipe inlet. 

on large-scale. As Re is slightly increased, the vorticity solution is considered to 
consist of two-term expansion, 

w(x,t) = 7 -b 21 /[ 7 ]^ 

The flow can still be regarded as laminar, and evidently the flow held has a large- 
scale structure because the space-time convolutions on 7 produce the flow scales 
comparable to the initial data. However, some increased degree of unsteadiness is 
detectable in the flow field due to the interaction among the eddies. If the Reynolds 
number is further increased, the instantaneous vorticity field is best described by 

a;(x, t) = 7 -I- 2 H [7 -I- 10 H [7 ]^, 


or possibly, by 


a;(x, t) = 7 -b 2 l/[ 7 ]2 -b 10 7 ]3 -b 62 7 

The water in the pipe now consists of multiple eddies of smaller scales which in¬ 
duce pronounced fluctuations in the flow. In practice, Reynolds did not observe the 
vorticity field but the velocity field. In view of the induction effect of the Green 
functions, the velocity held must show a space-time intermittent character due to 
the localised vorticity patches by the interaction of the interwoven 7 and the wall 
layer (cf. ipb)- Reynolds referred the regions of the velocity induced by the vorticity 
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r Experiments (Durst & Uensal 2006) 



Figure 5. In this experimental investigation for natural transition in a pipe of diameter ISmni 
and a length x = 10m, constant mass flows via a feedback control were maintained. To a 
very good approximation, a Hagen-Poiseuille profile was achieved at every Reynolds num¬ 
ber i?e<10000. Transition was detected at 7?etr~13000. At i?e=2030, 5000,10030,13000, 
the values of vt/cP' were 1.313,0.533,0.266,0.205 respectively. All the profiles were mea¬ 
sured at the pipe exit. The manner in which the flow was initiated suggested that the 
profile at the inlet ought to be either a constant or a boundary layer. It is seen that, at 
i?e=10030, the measured flow was an evolving profile. As Re was increased to Retr, the 
transitional profile was certainly non-self-similar. The test data at i?e=15020 and 16980 
showed that the process of losing laminarity initiated and finalised at locations upstream 
of the exit. In fact, the local flows over transitional process need not to be any partic¬ 
ular type as long as its progenitor at start is sizeable such that the non-linearity in the 
governing equations is evoked. 

concentrations as flashes (now commonly known as turbulent puffs and slugs). The 
whole flow field evolves into a structure which looks like a random motion. The 
Green functions is global and takes into account the downstream transitional or 
weakly turbulent flow state which in turn randomises the upstream localised vor- 
ticity slugs. At a higher i?e, the number of smaller eddies grows rapidly so that the 
flow becomes 

a;(x,t) = 7 -b 21 /[ 7]2 -f 101/[7]3-h621/[7]^-b430y[7]'^-f 3194y[7]®. 

The vorticity held must appear to be chaotic and irregular due to the strong inter¬ 
actions among the vorticity eddies. As Re is increased further, numerous gradated 
eddies emanate out of the non-linearity and viscous dissipation intensihes. The ran- 
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domness mechanism becomes marked and effective. The whole flow soon becomes 
chaotic as well as random as it contains a multitude of eddies in distinct scales. 
Beyond this stage, the flow field develops into full turbulence. 

Turbulence transition in pipe flow is an intrinsic, three-dimensional, non-linear 
process during the evolution. As the Reynolds number increases, the orderly stream¬ 
lined structure of initial laminarity is successively modified by the proliferated eddies 
of progressively smaller size until the entire structure is dissolved into a composi¬ 
tion of fluctuating shears. Depending on given initial data, there exists a Re-number 
range in which the non-linearity first becomes conspicuous. The process of transition 
can be an abrupt event or develops intermittently over a large section of space-time. 

(6) Turbulent puffs and slugs 

Local turbulent puffs are highly interactive vorticity due to the coupling of 
the r and 9 components as given in the governing equations (j3.15|) . The resolvent 
kernel rooted in (15.911 or more appropriate (16.1011 accumulates the components into 
local vorticity blobs. The convolutions ought to be stronger for the vorticity in 
wall layer. The concentrated vortices is advected downstream by the mean stream- 
wise velocity w so that they become elongated and spread out, forming localised 
vorticity as puffs. The process of proliferation by the non-linearity tends to weaken 
the local vorticity structures and enhances the local vorticity scales. Turbulent slugs 
are amalgamated puffs due to the scale multiplications. The puffs and slugs exist 
locally over part of pipe flow, it is reasonable to expect that the resulting motion 
incurs intermittent characters due to the induction of the vorticity in the whole 
pipe. 

Turbulent puffs refer to localised patches of disordered flow, typically 10—20 
diameter long. They are normally found at 1600<i?e<2500 in good experimental 
set-ups and advect downstream at 90—95% of the pipe mean velocity (see, for in¬ 
stance, Wygnanski & Champagne 1973). They are often called “equilibrium puffs” 
because they have structured flow features which do not seem to vary among dif¬ 
ferent experiments. In most experiments, the deductions of their lifetime or the 
travelling speed must be inferred from the ensemble-averages of the measurements. 
Our theory suggests that in a typical puff the local vorticity consists of flow scales 
in the order of 

10 ) 7 ]^ to 430 ( 7 ]^ (turbulent puffs). 

The mean sizes of the vortices ought to be moderate and are well-above the dis¬ 
sipative scales. Thus the turbulent puffs are not random in general, as compared 
to turbulence, but their dynamic behaviours do require statistical descriptions as 
the number of vorticity scales involved is already rather large. It is essential to 
notice that these puffs are not strong and abundant enough to trigger transition 
as they represent the transitional flow structures between the streamlined laminar 
flow and weakly turbulence. As the initial data become stronger or the Reynolds 
number becomes larger, more vortices of smaller scales will be created out of the 
non-linearity. As indicated in the re-laminarisation of the puffs by reducing the 
Reynolds number (Peixinho & Mullin 2006), the puff structures in the subsequent 
development remain in an identifiable fashion and they do not become chaotic. In 
particular, the test results suggest that the decay of the puff follows an exponential 
rate in space and in time (a typical smoothing effect due to heat kernel). The flow 
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visualisation highlights the existence of two longitudinal vortices preceding the ap¬ 
pearance of the puff flow at Reynolds number i?e=1750. It can be interpreted from 
many experimental data that, at the Reynolds number i?e<1750, the pipe flow has 
large-scale non-chaotic structures (§2 of Mullin 2011). The two terms solution of 
(EB) precisely describes the flow composition in the early evolutional phase. 

Turbulent slugs is related to the next flow scales in the vorticity hierarchy of 
the non-linearity. They can be described as vorticity patches with the scales from 
the vorticity powers 

430 ) 7 ]® to 3194 ) 7 ]® (turbulent slugs). 

The slugs are more space-filling as demonstrated in flow visualisations. The fluid 
motions in slugs show strong agitations and contain localised hysteresis due partly 
to the downstream flow effects. However, at Reynolds number i?e>3000, the flow 
state of initialisation turbulence is characterised by the presence of slugs. The sta¬ 
tistical properties of slugs ought to have similarities typically found in turbulence. 
However, the flow inside slugs is unlikely to resemble full turbulence. Unfortunately, 
current experimental techniques may not be able to uncover every detail of slug 
flow, given the presence of the numerous vortices. Hence use of statistical means 
for transition study is inevitable. 

(c) The intrinsic transition 

The ultimate aim of transition investigation is to identify the process or the 
mechanism which instigates the numerous vortices of various scales, apparently out 
of orderly streamlined progenitor flow. In reality, experimental evidence shows a 
great variety in the actual processes, depending on the types of the flow in ques¬ 
tion (free shear flows, boundary layers, confined flows in pipes or closed channels). 
The geometry of circular pipe is simple for experiments though we have seen that 
the analytic solutions are much more complicated compared to free shear flows in 
unbounded space. It is known that pipe flow development is experiment-dependent 
because of varied test conditions. Transition phenomena in different test environ¬ 
ments can be compared only if the disturbance levels are kept below some threshold 
values (see below §[5]). Should disturbances be significant and potentially aberrant, 
flow gradations to turbulence, as an initial-boundary problem of the equations of 
motion, may take divergent paths and case-dependent routes. We call the laminar- 
turbulent transition which is free of influence of disturbances the intrinsic transition. 

Proliferating vortices does not necessarily occur over all cross-sectional area 
at any hxed stream-wise z location. The intermittency is a natural phenomenon 
during flow evolution and transition. Whether local intermittent characters will 
persist into the subsequent turbulence depend on the initial data and the vicinity 
flow conditions. 

In practice, velocity profiles other than the self-similar distribution can also be 
generated and maintained. In figure [5] we plot and compare selected test data of 
Pfenninger (1961). It is comprehensible that delayed transition up to Re^lO® was 
achieved. These tests ought to be assessed as dedicated flow controls to postpone 
turbulence transition in the same way as suction is applied to laminar boundary 
layers in high Reynolds number. In theory, the test Ve must be so powerful that 
only the hrst two terms in the vorticity series are adequate to describe the initial 
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Experiments (Pfenninger 1961^ 


0.9 - 


R^ = 0.00156 
= 0.00441 
R^ = 0.00623 
R^ = 0.03094 

Cl 

— Parabolic profile 


Figure 6. Measured velocity profiles in accelerating flows in favourable pressure gradients 
at various Rd = 4:{x* jd) jRe, where pipe diameter d — 50.8mm, and x* is the length of an 
equivalent straight pipe with zero boundary layer thickness at the entry. For the values 
of Rd, the axial distances from the entry are 1.37m, 3.045m, 5.94m and 11.88m respectively. 
The velocity profiles substantially deviate from the parabolic profile over the whole pipe 
length. The main experiments were conducted in accelerating flows, likely involving strong 
swirls. The energy of the flows in the favourable pressure gradients was much higher than 
the parabolic profile in constant pressure gradient at the same Reynolds number. 


evolution over the complete pipe length. On the other hand, there does not ex¬ 
ist experimental evidence to indicate that fully developed Hagen-Poiseuille flow of 
parabolic profile remains laminar at arbitrary Reynolds number. 

On the basis of the outcomes from credible constant-mass-flow-rate experiments 
(see, for example, Darbyshire & Mullin 1995; Durst & Uensal 2006), a number of 
conclusions are compatible with the regularity theory of the Navier-Stokes equa¬ 
tions. An array of test techniques, including suction, jets, ring blockage and iris 
diaphragm, has been utilised to instigate finite perturbations into flows at high 
Reynolds number. The established results showed that transition is relatively in¬ 
sensitive to the means of forcing which are the representation of initial data. In par¬ 
ticular, good repeatability in transition measurements has been confirmed in both 
natural and forced scenarios. One may opt the term in every component of (ED 
as the unknown in view of u=A the equations then become a set of quasi-linear 
differential equations, instead of non-linear ones. Such a view may be instructive 
in understanding the quantitative repeatability of the test results, that have been 
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collected at the Reynolds numbers preceding the full-blown random effect due to 
viscous dissipation. The differential operators, and dt—vA' in (I3.15L are 

the stereotype of continuum diffusion process for finite v>Q. The flow evolution 
governed by the Navier-Stokes equations does not behave like certain lower-order 
dynamical systems which are sensitive to initial conditions. 

A theoretical development of pipe flow transition is to identify particular numer¬ 
ical solutions known as travelling waves. They are identified as unstable solutions 
of the governing equations and are stream-wise vortices propagating with constant 
phase speed. Specifically, the wavy topology has been considered to originate from 
saddle-node bifurcations of the stability boundaries in the Reynolds number range 
770<i?e<1250 (Faisst & Eckhardt 2003; Wedin & Kerswell 2004). However, the 
Navier-Stokes equations are deterministic and globally well-posed; exact solutions 
of the governing equations do not bifurcate in space or in time. 

8. Experimental evidence for disturbance-free transition 

One of the key outcomes from the present theoretical analysis is that the laminar- 
turbulent transition is a direct consequence of the non-linear term (u.V)u in the 
equations of motion. All the integral power terms in (j7.1ll are emanated from the 
non-linearity. Thus our analysis asserts that the intrinsic transition occurs even in 
the complete absence of disturbances, infinitesimal or finite. Since fluid motions 
are an initial-boundary value problem of the Navier-Stokes equations, initial and 
boundary conditions make allowance for any disturbances. 


Microscopic fluctuation and turbulence 

One may contemplate that microscopic fluctuations are the internal disturbances 
that play a role in turbulence transition. However, molecular agitations are incessant 
in all types of fluid motions; the fluctuations exist in laminar, transitional as well 
as turbulent flows. In fact, their existence was first discovered in fluids at rest 
(Brownian motion). Turbulent flows in finite-energy motion are dynamic in nature 
and hence do not persist indefinitely. Experiments demonstrate that turbulence in 
pipe flow undergoes reverse transition during flow decay when a large share of the 
energy has been dissipated. In most laboratory tests under standard temperature 
and pressure, any temperature increase from the viscous dissipation of flow energy 
is not critical enough to alter the mean molecular speed according to kinetic theory 
of gases (Maxwell 1867; Boltzmann 1906). Eor incompressible flows, the internal 
energy of fluid’s molecules has a primarily identical Boltzmann distribution in all the 
three flow stages. If the microscopic fluctuations are responsible for the transition 
from laminar to turbulent state, in the absence of any macroscopic disturbances, 
are they the cause for the process of reverse transition? 

Nevertheless, the microscopic fluctuations are identified as the origin of the 
random characters in turbulence in incompressible flows (Lam 2013). Viscous dissi¬ 
pation on the copious eddies of small scales expedites entropy production in order 
to restore the flows to the equilibrium of shear-free state. The diffusion of individual 
small-scale eddies is a stochastic process in space and in time as fluid’s thermal en¬ 
ergy unceasingly fluctuates. Fluid motions on the continuum inhere the microscopic 
randomness through the diffusive dissipation process. 
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Environmental disturbances 

These are the external disturbances, such as tree-stream turbulence, acoustic 
noise, mechanical vibrations of test apparatus. They can be specified in the initial 
data (I3.11|l or (|3.19L and in the boundary value (I3.7|) . In principle, their effect on 
transition has been fully taken into account in the equations of motion. Specifi¬ 
cally, any background disturbance at the far end of the pipe is fully enumerated 
in the Green function (ED or (16.31) which encloses the entire pipe interior. The 
Fourier-Bessel expansion in the function explicates the flow field as a space-bounded 
“wave-like” structure with diminishing magnitude as vt^oo. The inherent analytic 
characters of the solution precisely explain many wavy patterns observed in flow 
visualisations (see, for instance, van Dyke 1982 and Samimy et al 2003). 

In the classic experiments of Schubauer & Skramstad (1947) for transition flows 
on a flat-plate in zero pressure gradient, it was demonstrated that the transition 
location is independent of the free-stream turbulence once the level dropped below 
a threshold about 0.08%. The turbulence level is defined as 

T^= (^J(^ + 1^ + 

expressed in percentage. It was speculated that the acoustic excitation in the tunnel 
might be responsible for the transition. Additional tests showed that the transition 
was delayed if environment noise was carefully regulated (Wells 1967). In particular, 
when the noise level was below 90db, the transition was basically insensitive to the 
ambient disturbances as long as Tu<0.1%. In figure[7]we summarise the transition 
location as a function of the free-stream turbulence with reduced acoustic excita¬ 
tion. The limiting case of zero disturbances can be obtained by extrapolation. In 
fact, these limits have been drawn in the early reviews, see, for example, figure 10 of 
Tani (1969) or figure 16(a) of Narasimha (1985). However the significance of these 
limits has not been realised due probably to lack of theoretical evidence at the time. 
Instead, the occurrence of the natural transition was attributed to residue or back¬ 
ground disturbances such as derivation from inlet flow, thermal convection due to 
temperature differences or imperfections in test apparatus. Nevertheless, all these 
perturbation irregularities can be minimised by means of control, and in principle 
they are best regarded as extraneous for fluid motions. 

We draw our attention to a puzzling fact in flat-plate experiments: in order 
to observe the Tollmien-Schlichting waves, artificial excitations of specific frequen¬ 
cies had to be introduced into the laminar boundary layer so that the resulting 
measurements would agree with the predictions by linear stability theory which 
is an attempt to solve the complete Navier-Stokes equations by approximation. In 
the complete absence of the simulated excitations, instability waves could not be 
detected at the neighbourhood of the critical Reynolds number as long as the ex¬ 
ternal environmental disturbances were well-kept below a minimum, even though 
turbulence transition was unavoidable at some downstream location (Schubauer & 
Skramstad 1947; Klebanoff et al 1962). An explanation of this apparently paradoxi¬ 
cal situation is that the Tollmien-Schlichting waves are usually too weak to uncover 
with certainty (see, for instance, §17.7 of Tritton 1988). The question is: how can 
we distinguish the man-made disturbances from the “genuine” instability waves? 
Evidently, if the flow has been disturbed, the measurement instrument ought to 
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Figure 7. Transition Reynolds number {Rextr) as a function of free-stream turbulence 
level Tu for flat-plate boundary layers. In each case, the intrinsic turbulence transition is 
found to be independent of the threshold value T* as shown by a small vertical line at 
0.1% and 0.08%. If the turbulence levels are substantially higher than T^, loss of laminar 
flow occurs prematurely but the actual process is mis-shaped and contorted by the test 
environment. Incoming flows with strong fluctuations can induce complete loss of laminar 
flow as the first two terms in our series solution (ESI contain convoluted vortex eddies of 
the free-stream turbulence. For the two perturbative environments commonly encountered 
in practice, disturbance-free transition is established at a finite value of Rextr (marked 
by a filled diamond) which is determined by extrapolating the test data to Tu = 0. 


have picked up the responses of the laminar flow to the forced excitations since, in 
any event, the unstable waves are ineffective to Impact. 

A comparable situation was realised in the leading-edge flow over a swept wing 
where the three-dimensional boundary layer was dominated by cross-flows (Deyhle 
& Bippes 1996). In particular, controlled excitations or artificial disturbances had 
to be generated in order to observe cross-flow instabilities as stipulated in linear 
stability theory. The test conclusion says that free-stream turbulence affects transi¬ 
tion indirectly as long as ru<0.2%. It is also found that the receptivity effects due 
to acoustic excitations and non-uniformities of mean flow are weak. The cross-flow 
instabilities are seen less vulnerable to free-steam turbulence but the trend is in 
accord with the flat-plate experiments. By analogy, the main results may be inter¬ 
preted as a confirmation that unstable cross-flow waves as predicted by the linear 
theory are unlikely observable in the low turbulence test environments. 
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To obtain a solution from the equations of motion, any spatio-temporal flow 
modification due to the presence of a disturbance during evolution requires revised 
initial data and a “new” solution of the governing equations must be sought. This is 
particular relevant to experiments where injections are used as a means of creating 
finite disturbances on pipe flow. Since the mass addition must be taken into account 
as V.u=/(a;)f^0, any deviation from the continuity hypothesis must be assessed lo¬ 
cally in a control volume analysis. The vorticity field is no longer solenoidal and a 
modification in the no-slip boundary condition has to be considered. As a result of 
the problem reformulation, there are extra terms in the equations of vorticity dy¬ 
namics which can be aggressive to reorganise the overall flow evolution. The forced 
phenomena observed in these related experimental work ought to be quantified in 
dedicated studies. 


9. Conclusion 

We have shown that the flow evolution in a circular pipe from given initial state 
can be described by solution of the Navier-Stokes equations. The general solution 
of the vorticity equation offers a theoretical explanation to the experimental ob¬ 
servations of Reynolds (1883). The non-linear effects specify a dynamic process of 
proliferating vorticity scales. We have also demonstrated that the vorticity produc¬ 
tion at pipe wall cannot be prescribed a priori, but can be determined by iteration. 
It is a matter of computations to map out the detailed flow structures. In essence, 
the vorticity solution not only explicates turbulence but also clarifies the space- 
time gradation of the flow evolution. Specifically, the laminar-turbulent transition 
is understood as an intrinsic part of the process should the initial vorticity or the 
Reynolds number be sufficiently large. It is inevitable for laminar, transitional and 
turbulent motions co-exist over a large part of the pipe so that the flow field is nor¬ 
mally intermittent in space and in time. In light of test evidence and our analytical 
results, we assert that the intrinsic turbulence transition is primarily instigated by 
the non-linearity. Disturbances, either infinitesimal or finite, misrepresent the tran¬ 
sition. On practical as well as theoretical grounds, it is genuinely conceivable to 
understand the transition dynamics from the primitive equations without recourse 
to modelling because exact vorticity solutions can be obtained. Consequently, the 
impact of the non-linear effects on the flow development can be appreciated with 
confidence. Given the growing intricacy in flow field prior to fully random turbulent 
state, experimental investigations on transition process are surely overwhelmed by 
aberrant particulars and hence have to, more or less, rely on statistical inference. 


28 May 2015 
f . lamllSyalioo. com 
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Appendix A. Zeros for Bessel functions Jn{o'n,k) and J'n{an,m) 

For every n, the zeros of the Bessel function = 0 may be arranged in 

ascending order 0 < an,i < <Jn ,2 < o'n ,3 < ■ ■ ■ < <Tn,k < • • • • The ascending order is 
also true for the zeros of J^(a„_i) = 0. As n —>■ oo, the zeros can be estimated by 
the asymptotic formula (Olver et al 2010) 

CT„,i ~ n + 1.8557571ni/3 + 1.03315/n^/3- 0.00397/n, 


d„,i ~ n + 0.808616571^/3+ 0.07249/ni/3- 0.05097/n. 

We list the first few zeros of Jn{(^n,i) = 0 and J'n{an,i) = 0 in the table below. 


n 


^n,l 

n 

Q^n,l 


0 

2.404826 

- 

0 

0.000000 

- 

1 

3.831706 

3.884937 

1 

1.841184 

1.830137 

2 

5.135622 

5.156134 

2 

3.054237 

3.050843 

3 

6.380162 

6.391488 

3 

4.201189 

4.199499 

4 

7.588342 

7.595682 

4 

5.317553 

5.316522 

5 

8.771484 

8.776696 

5 

6.415616 

6.414913 

6 

9.936110 

9.940037 

6 

7.501266 

7.500752 

7 

11.086370 

11.089456 

7 

8.577836 

8.577441 

8 

12.225092 

12.227593 

8 

9.647422 

9.647107 

9 

13.354300 

13.356376 

9 

10.711434 

10.711176 

10 

14.475501 

14.477256 

10 

11.770877 

11.770661 

11 

15.589848 

15.591356 

11 

12.826491 

12.826308 

12 

16.698250 

16.699562 

12 

13.878843 

13.878685 

13 

17.801435 

17.802589 

13 

14.928374 

14.928237 

14 

18.899998 

18.901023 

14 

15.975439 

15.975317 

15 

19.994431 

19.995348 

15 

17.020323 

17.020215 

16 

21.085146 

21.085973 

16 

18.063265 

18.063168 

17 

22.172495 

22.173244 

17 

19.104462 

19.104375 

18 

23.256776 

23.257459 

18 

20.144083 

20.144003 

19 

24.338250 

24.338876 

19 

21.182270 

21.182197 

20 

25.417141 

25.417717 

20 

22.219146 

22.219080 

21 

26.493647 

26.494180 

21 

23.254821 

23.254759 

22 

27.567944 

27.568438 

22 

24.289386 

24.289329 

23 

28.640185 

28.640644 

23 

25.322924 

25.322872 

24 

29.710509 

29.710937 

24 

26.355510 

26.355461 

25 

30.779039 

30.779440 

25 

27.387207 

27.387161 

26 

31.845887 

31.846264 

26 

28.418075 

28.418032 

27 

32.911154 

32.911508 

27 

29.448165 

29.448125 

28 

33.974930 

33.975264 

28 

30.477526 

30.477488 

29 

35.037299 

35.037614 

29 

31.506199 

31.506163 

30 

36.098337 

36.098635 

30 

32.534224 

32.534190 
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